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' the dust epoch, without a fine tunning because the contribution of the scalar field in this 

theory is equivalent to a stiff fluid (as can be seen from the energy-momentum tensor for 
the scalar field), that evolves in a different way as the dust component. To have similar 
contributions of the scalar component and the dust component implies that their past values 
were fine tunned. So, we reinterpreting this null result as an indication that dark matter 
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effects in the cosmic microwave bound radiation, and than this formalism yield to this epoch 
as primigenius results. We do the mention that this formalism was used recently in the 
so called K-essence theory applied to dark energy problem, in place to the dark matter 
problem. Also, we include a quantization procedure of the theory which can be simplified by 
reinterpreting the theory in the Einstein frame, where the scalar field can be interpreted as 
part of the matter content of the theory, and exact solutions to the Wheeler-DeWitt equation 
are found, employing the Bianchi Class A cosmological models. 
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I. INTRODUCTION 

Several observations suggest that in galaxies and galaxy clusters there is an important quantity 
of matter that is not interacting electromagnetically, but only through gravitation. This is the well 
known dark matter problem. Several solutions have been consider for this problem, modifying the 
gravitational theory or introducingnew forms of matter and inter accions. To adress the dark matter 
problem Saez and Ballester (SB) [l| formulated a scalar-tensor theory of gravitation in which the 
metric is coupled with a dimensionless scalar field. In a recent analysis using the standard scalar 
field cosmological models , contrary to claims in the specialized literature, it is shown that 
the SB theory cannot provide a realistic solution to the dark matter problem of Cosmology for 
the dust epoch, because the contribution of the scalar field is equivalent to stiff matter. We can 
reinterpret this result in a sense that the galaxy halo was formed during this primigenius epoch and 
its evolution until the dust era using the standard scalar field cosmological theory. In this theory 
the strength of the coupling between gravity and the scalar field is determined by an arbitrary 
coupling constant lv. This constant uj can be used to have a lorenzian (-1,1,1,1) or seudo-lorenzian 
(-1,-1,1,1) signature when we build the Wheeler-DeWitt equation. The values for this constant, in 
the classical regime, are dictated by the condition to have real functions. Other problem inherent 
to this theory is that not exist how build the invariants with this field as in the case to scalar 
curvature. So, was necessary to reinterpret the formalism where this field is considered as matter 
content in the theory in the Einstein frame. 

On the other hand, this approach is classified with another name, by instant, Armendariz- 
Picon et al, called this formalism as K-essence ^j, as a dynamical solution for explaining naturally 
why the universe has entered an epoch of accelerated expansion at a late stage of its evolution. 
Instead, K-essence is based on the idea of a dynamical attractor solution which causes it to act as a 
cosmological constant only at the onset of matter domination. Consequently, K-essence overtakes 
the matter density and induces cosmic acceleration at about the present epoch. Usually K-essence 
models are restricted to the Lagrangian density of the form 

:|4 



d^x^m {v<pr. (1) 

n 

One of the motivations to consider this type of Lagrangian originates from string theory [5]. For 
more details for K-essence applied to dark energy, you can see in [S] and reference therein. Many 
works in SB formalism in the classical regime have been done, where the Einstein field equation 
is solved in a direct way, using a particular ansatz for the main scalar factor of the universe 



B y, Q, [n 



l|, yet a study of the anisotropy behaviour trough the form introduced in the hne 



element has been conected 



a, 
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On another front, the quantization program of this theory has not been constructed. The main 
comphcation can be traced to the lack of an ADM type formalism. We can transform this theory 
to conventional one where the dimensionless scalar field is obtained from energy-momentum tensor 
as an exotic matter contribution, and in this sense we can use this formalism for the quantization 
program, where the ADM formalism is well known 21 1. 



In this work, we use this formulation to obtain classical and quantum exact solutions to 
anisotropic Bianchi Class A cosmological models with stiff matter. The first step is to write SB 
formalism in the usual manner, that is, we calculate the corresponding energy-momentum tensor to 
the scalar field and give the equivalent Lagrangian density. Next, we proceed to obtain the corre- 
sponding canonical Lagrangian Ccan to Bianchi Class A cosmological models through the Legendre 
transformation, we calculate the classical Hamiltonian from which we find the Wheeler-DeWitt 
(WDW) equation of the corresponding cosmological model under study. We employ in this work 
the Misner parametrization due that a natural way appear the anisotropy parameters to the scale 
factors. 

The simpler generalization to Lagrangian density for the SB theory [l| with the cosmological 
term, is 

^geo = (R - 2A - F((^),/.,^</.'^) , (2) 

where (p''^ = g'^^cp^a, R the scalar curvature, F{(j)) a dimensionless function of the scalar field.. In 
classical field theory with scalar field, this formalism corresponds to null potencial in the field (j), 
but the kinetic term is exotic by the factor F{(j)). 

From the Lagrangian ([2]) we can build the complete action 

1= / V^iCg^o + C^^t)d\ (3) 
Jt, 

where Cmat is the matter Lagrangian, g is the determinant of metric tensor. The field equations 
for this theory are 

Ga(3 + ga/3A - F{<p) (^(p^a<P,f} - ^ga/3'/',7<^''^^ = -SttGTo^, (4a) 

2F(0)<A'f, + = 0' (4b) 

where G is the gravitational constant and as usual the semicolon means a covariant derivative. 
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The equation ()4bp take the following form for all cosmological Bianchi Class A models, assuming 
that the scalar field is only time dependent ( here / = ^ = -j^) 

3n'ct>'F + ct>"F + l^^'' = 
2 ag) 

which can be put in quadrature form as 

\f^'' = Foe-6^, (5) 

this equation is seen as corresponding to a stiff matter content contribution. 

The same set of equations (|4al4bp is obtained if we consider the scalar field (j) as part of the 
matter budget, i.e. say C(f, = —¥{(f))g°'^(j)^a4>^i3 with your corresponding energy-momentum tensor 

= F(0) (^^,a<P,p - \gaP<Pa^'^^ ■ (6) 

which is conserved directly considering a stiff matter era in a barotropic scalar fluid (see appendix 
section 8). In this new line of reasoning, action ([3]) can be rewritten as a geometrical part (Hilbert- 
Einstein with A) and matter content (usual matter plus a term that corresponds to the exotic 
scalar field component of SB theory). 

In this way, we write the action ([3]) in the usual form 

I = £ ^/=g (R - 2A + £^at + C-^) d^x, (7) 

and consequently, the classical equivalence between the two theories. We can infer that this corre- 
spondence also is satisfied in the quantum regime, so we can use this structure for the quantization 
program, where the ADM formalism is well known for different classes of matter [21]. Using this 
action we obtain the Hamiltonian for SB. We find that the WDW equation is solved when we 
choose one ansatz similar to this employed in the Bohmian formalism of quantum mechanics and 
the gravitational part in the solutions are the same that these found in the literature, years ago 
22|. 

This work is arranged as follow. In section 3 we present the method used, employing the FRW 
cosmological model with barotropic perfect fluid and cosmological constant. In section 4 we con- 
struct the Lagrangian and Hamiltonian densities for the anisotropic Bianchi Class A cosmological 
model. In section 5 the classical solutions using the Jacobi formalism are found. Here we present 
partial results in the solutions for some Bianchi's cosmological models. Classical solution to Bianchi 
I is complete in any gauge, but the Bianchi II and VIh=-i, the solutions are found in particular 
gauge. Other Biachi's, only the master equation are presented. In Section 6 the complete cuantiza- 
tion scheme is presented, obtaining the corresponding Wheeler-DeWitt equation and its solutions 
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are presented 
and Shepley, 



in unified way using the classification scheme of Ellis and MacCallum 33] and Ryan 



33| 



II. THE METHOD 

Let us start with the line element for a homogeneous and isotropic FRW universe 



-N\t)dt' + a'{t) 



1 — Kr^ 



(8) 



where a{t) is the scale factor, N{t) is the lapse function, and k is the curvature constant that can 
take the values 0, 1 and —1, for flat, closed and open universe, respectively. The total Lagrangian 
density then reads 



N 



6KNa + 



N 



+ 167rGiVa> - 2Na^A , 



(9) 



where p is the energy density of matter, we will assume that it complies with a barotropic equation 
of state of the form p = 7p, where 7 is a constant. The matter content is assumed as a perfect 
fluid Tfj_u = (/9 + p)Ufj,Uu + g^ii^p where is the fluid four-velocity satisfying u^u'^ = —1 . Taking 
the covariant derivative we obtain the relation 

snp + Stlp + p = 0, 

whose solution becomes 



p^e 



-3n(i+7) 



(10) 



where p-y is an integration constant. 

From the canonical form of the Lagrangian density ([9]), and the solution for the barotropic 
fluid equation of motion, we find the Hamiltonian density for this theory, where the momenta are 
defined in the usual way Ugi = where q^ = (a, </)) are the field coordinates for this system. 



dC 

da 



12aa 

IT' 



dC 2Fa^<j) 

= — - -- 



12a ' 



2Fa^ ' 



(11) 



so, the Hamiltonian density become 



n 



24 



a^Ul + 



F{<P) 



Ul + 144Ka^ + 48a^A - 3847rGp^a 



3(1-7) 



(12) 
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IS 



dS ______ 

Using the transformation liq = the Einstein-Hamilton- Jacobi (EH J) associated to Eq. 



The EHJ equation can be further separated in the equations 



6 / dS^, \ 2 



2 

= /^', (14) 



(^) + 48a'A - 3847rGp^a^^^-^^ = -^^ , (15) 

where is a separation constant. With the help of Eqs. (|lip . we can obtain the solution up to 
quadratures of Eqs. and (fTSj) . 

I y^)dct> = ^1 a-Hr)dT, (16a) 
Ar = / , (16b) 
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with = Eq. (I16ap readily indicates that 

F((^)(/.2 = 6z^2^-6(^)^ (17) 

Also, this equation could be obtained by solving equation ()4bp . Moreover, the matter contribution 
of the SB scalar field to the r.h.s. of the Einstein equations would be 

= ^F(0),/.2 oc . (18) 



this energy density of a scalar field has the range of scaling behaviors [23l. I24i| . is say, scales exactly 
as a power of the scale factor like, oc a""*, when the dominant component has an energy density 
which scales as similar way. So, the contribution of the scalar field is the same as that of stiff 
matter with a barotropic equation of state 7 = 1. This is an interesting result, since the original 
SB theory was thought of as a way to solve the missing matter problem now generically called the 
dark matter problem. To solve the latter, one needs a fluid behaving as dust with 7 = 0, it is 
surprising that such a general result remains unnoticed until now in the literature about SB. This 
is an instance of the results of the analysis of the energy momentum tensor of a scalar field by 



Marden [251] for General Relativity with scalar matter and by Pimentel [2a] for the general scalar 
tensor theory. In both works a free scalar field is equivalent to a stiff matter fluid. 

Furthermore, having identified the general evolution of the scalar field with that of a stiff fluid 
means that the Eq. (flBb]) can be integrated separately without a complete solution for the scalar 
field. In ^ appear a compilation of exact solutions in the case of the original SB theory to FRW 
cosmological model and in ^] were presented the classical and quantum solution to Bianchi type I. 
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III. THE MASTER HAMILTONIAN TO BIANCHI CLASS A COSMOLOGICAL 

MODELS 

Let us recall here the canonical formulation in the ADM formalism of the diagonal Bianchi 
Class A cosmological models. The metric has the form 

= -dt^ + e^m {e^f'^^ju'cJ, (19) 

where /3ij(t) is a 3x3 diagonal matrix, /3ij = diag(/3+ + \/3/5-,/3+ — \/3/?-, —2(3^), fi{t) is a scalar 
and are one-forms that characterize each cosmological Bianchi type model, and that obey du^ = 
^Cjjja;j A cu^, Clj^^ the structure constants of the corresponding invariance group, these are included 
in table 1. 



Bianchi type 


1-forms ui^ 


I 


uj^ = dx^, Lo'^ = dx^, Lo''' = dx^ 


II 


uj^ = dx^ - xMx'"^, = dx-^, a;^ = dx^ 


VIh=-i 


ui^ = e'^'^dx^, cj^ = e'^Mx^, co^ = dx^ 


VIIo 


uj^ = dx^ + dx^, a;^ = — dx^ + dx^, uj^ = dx^ 


VIII 


CO^ = dxl + [1 + (xl)2]dx2 + [xl - x2 - (xl)2x2]dx3, 

oj^ = 2xMx2 + (1 - 2xix2)dx=^, 

oj^ = dx^ + [-1 + {x'^fjdx'^ + [x^ +x^ - {x'^fx'^]dx^ 


IX 


uj^ = — sin(x^)dx"'^ + sin(x^) cos(x^)dx^, 

= cos(x^)dx^ + sin(x^) sin(x^)dx^, = cos(x^)dx^ + dx^ 



Table 1. one- forms for the Bianchi Class A models. 



We use the Bianchi type IX cosmological model as toy model to apply method discussed in the 
previous section. The total Lagrangian density then reads 

i2 a2 a2 



Pi 

N N N 



N ^ 



+ 167rGNp - 2NA 



_ ^g-2/3++2V3^_ g-2^+-2^/3/3_ _|_ ^4^+^ | 



(20) 
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making the calculation of momenta in the usual way, IlqM = where q'^ = (Q, /J+j /?_, i 



N ' 12 



and introducing into the Lagrangian density, we obtain the canonical Lagrangian as 

with the Hamiltonian density 

p-3n / fi \ 

^± = ^ (-n^ - + + n^. + u(i7, /3±) + Ci J , (21) 

where the gravitational potential becomes, 
U(0, /3±) = 12e^^ + ^^^^-^Vs/S- + ^4/3+ _ 2 ^^^m + e2'^+-2V3/3_ ^ g-2/3++2V3/3_ 1^ ^ 

with Ci = 3847rGpi corresponding to stiff matter epoch, 7 = 1. 

The equation (I2ip can be considered as a master equation for all Bianchi Class A cosmological 
model in the stiff epoch in the Saez-Ballester theory, with \J{Q,f3±) is the potential term of the 
cosmological model under consideration, that can read it to table II. 



IV. CLASSICAL SCHEME 

In this section, we present the classical solutions to all Bianchi Class A cosmological models 
using the appropriate set of variables, 

^3 = 1^-2/3+. (22) 
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Bianchi type 


Hamiltonian density % 


I 


— 

C 

24 




fnl + + n?. - 48Ae6f^ + 3847rGp^e-3(7-i)f^l 


II 


24 


- f n2 + n2 + n?. - 48Ae6^ + 3847rGp^e-3(T-i)^ 
M2e^^e^''++^^'^-l 


VI_i 


24 


-n2 - f n2 + n2_ + - 48Ae6f^ + 3847rGp^e-3(7-i)i^ 
f48c4"c4^-*+] 


VIIo 


— 3 Si 

e 

24 


-n2 - 


+ n2_ + n?. - 48Ae"^ + 3847r( 

'g4/3++4\/3/3_ _ g4/3+ ^_ g4/3+-4v^/3_ J 




VIII 


24 


M2e4^2 
-2 {e4^ 


f n| + + n?^ - 48Ac6^^ + 3847rC 

'g4/3++4\/3/3_ _|_ g4/3+-4\/3/3_ _|_ g-8/3+ 
t _ g-2/3+-2\/3/3_ _ g-2/3++2V3/3_ |A 


5p^e-3(7-i)^^ 


IX 


24 


-n^ - f + n2_ + n?. - a^Kc^^ + 3847: 

_2 |e4/3+ + e2/3+-2\/3/3- + g-2/3++2V3/3_ jj" 


Gp^e-3(^-i)^ 



Table II. Hamiltonian density for the Bianchi Glass A models. 



A. Bieinchi I 



For building one master equation for all Bianchi Class A models, we begin with the simplest 
model give by the Bianchi I, and give the general treatment. The corresponding Lagrangian for 
this cosmological model is written as 

"2/3i/32 , 2/ji/33 , 2/52/33 , F(<^)</<2 



Ci = e' 



+ 



+ 



+ 



N N N N 

the momenta associated to the variables {^i , 0) are 



+ 16N7rGp^e-(^+^)('^i+^2+/33) _ 2NA 



^2 = ^e-(^i+^2+/'3)(n, + n3-n2), 



n2 = ^(/jl + /j3)e^^+^^+^^ 
n3 = ^01+ /32)e^^+^^+^^ /33 = J e-(/^^+/^^+'33)(ni + n2 - Ha), 



N 



iip-(/3i+/32+/33)n, 
2F ^ 



so, the Hamiltonian is 



1 



-uj - ni - ni + + 2nin2 + 20103 + 20203 



+16Ae2(/5i+/32+/33) _ i287rGp^e(i-^)('3i+^^+^3) 



(23) 



(24) 



(25) 
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n; = 


-4Ae^i+^2+/33 + i6^G(l 


- 1)p^e 


~7(/3i+/32+/33) 


n'2 = 


_4Ae/5i+/32+/33 + i67rG(l 


- l)p-l^~ 


-7(/3i+/32+/33) 


4 = 


-4Ae^i+^2+/33 + i6^G(l 


- l)p-i^~ 


-7(/3i+/32+/33) 


n; = 


1 TT' 

-^p-(/3i+/32+/33) * n2 

4 






/3l = 


^e-(/3i+&+/33) [_ni + n2 


+ n3], 




02 = 


ie-(/3i+&+/33) [_n2 + ni 


+ n3], 




= 


le-(/3i+&+/33) [_n3 + Hi 


+ n2], 




^ = 


J_p-(/3i+/32+/33)n^ 

2F 







using the hamilton equation, where / = ^ = we have 

(26) 
(27) 
(28) 
(29) 

(30) 
(31) 
(32) 
(33) 

equations (I26I27I28D implies 

Hi = n2 + ki = n3 + k2. (34) 

Also, the differential equation for field (j) can be reduced to quadratures when we use equations 
(1291) and ([Ml), as 

1f((^)(/.'2 =,/.oe-2(^i+/^2+/33)^ ^ y/fX^d^ = V^e-^'^i+'^^+^^^dr, (35) 

which correspond to equation ([5]) obtained in direct way from the original Einstein field equation. 
The corresponding classical solutions for the field (p for this cosmological model can be seen in ref. 

a. 

Using this result and the equation for the field (p given in (j24p we can find that 2^ = 1600- 
From the hamilton equation for the momenta Hi can be written for the two equations of state 
7 = ±1, introducing the generic parameter 

{-4A 7 = 1 

(36) 
-4A + 327rGpi 7 = -1 

as n'^ = XQP^+P^~^P'-i ^ then re-introducing into the Hamiltonian equation (|25p we find one differential 
equation for the momenta Hi as 

^n;2 + 2n2 - - kg = 0, (37) 

A 



where the corresponding constants are 



, k? + k^ - 16(/)o, 7 = -1 

Ac = 2(ki+k2), k3=<; ^ ' (38) 

kf + k| - 1600 + 1287rGpi , 7 = 1 
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and whose solution is 



K + 12k3 . 
Hi = - ± sm 



3A 



At 



(39) 



On the other hand, using this result in the sum of equation ()30|31|32p . we obtain that 

/3A 



/3i + /32 + /33 = Ln 



solution previously found in ref, 



.B 



a 



a2A 



COS 



At 



a = + 12k3 



(40) 



using the Hamilton-Jacobi approach. 



B. Bianchi's Class A cosmological models 



The corresponding Lagrangian for these cosmological model are written using the Lagrangian 
to Bianchi I, as 



Cii = A + Ne^i+^2+/33 
/:vii,.o = + Ne^i+^2+^« 



1 



,2(/3i-/32-/33) 



2 

2e-2fe 
1 



1 



2(/3i-/32-/33) I ie2(-/3i+/32-/33) _ g-Sfe 

2" 2 

g2(/3i-/32-/33) _^ g2(-/3i+/32-/33) _^ g2(-/3i-/32+/33) 



(41) 
(42) 
(43) 

(44) 
(45) 



the momenta associated to the variables are the same as in equation (j64p . so, the generic 

Hamiltonian is 



2(-e-2/3i+e-2/52+e 



-2/33 



Ax = £1 + ^e'^^+z'^+Z^^ 



2 ( e-^/^i + e-^'^2 + g 



g2(/3i-/32-/33) _^ g2(-/3i+/32-/33) _^ g2{-/3i-/32+/33) 



2/32 , ^-2/33 



nA = ni- ie-(^i+^^+^3) [Ua(/?i,/32,/33)] , 



(46) 



where the potential term Ua(/3i, /32; /^s) is given in table III, where A corresponds to particular 
Bianchi Class A models (I, II, VIh=-i,VIIh=o,VIII,IX). If we choose the particular gauge to the 
lapse function N = eil^^+l^^+l^s) ^ ^j^jg equation (j46p is much simpler, 



^A = ^i- 7t[Ua(/3i,/32,/33)] 



(47) 



where Hi is as in equation (j25p but without the factor e (/^i+z's+fe) 
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Bianchi type 


Potential Ua(/3i, /32, /Ss) 


I 





II 


e4/3i 


VIh=-i 




VIIh=o 




VIII 




IX 





Table III. Potential \Ja{Pi, (^2, f^a) for the Bianchi Class A Models. 



The Hamilton equations, for all Bianchi Class A cosmological models are as follows 

n; = -4Ae'^i+'^2+/33 + i67rG(l-7)/>-,e-^('^i+^2+'^3) 

n'2 = -4Ae'^i+'^2+^3 + 167rG(l-7)/9^e-^('^i+^2+'^«) 
d 











-(/31+/32+/33) [_n, + n2 


+ n3 


= r 


-(/^i+fe+A) [_n2 + Hi 


+ n3 


« = r 


-(/31+/32+/33) [_n3 + m 


+ n2 


0' = ^ 

^ 2F 







(48) 



+^ (^^e-(''i+^^+'^3) [UA(/3i,/32,/33)]j , (49) 
n'3 = -AKe^^+P^+^i + l67rG(l - -^)p^e-^^^^+^^+^^^ 

^Ws (r"^"'^"'"'"'^ [Ua(/3i,/32,/33)]) , (50) 

(51) 
(52) 
(53) 
(54) 
(55) 

In this cosmological models, it is remarkable that the equation for the field (f) (j35p is mantained 
for all Bianchi Class A models, and in particular, when we use the gauge N = e^i-+l^2+l33 ^ ^j^jg 
solutions for this field are independent of the cosmological models. 

C. Classical solution in the gauge N = e^i+'^^+^^j A = and 7 = 1 



With these initial choices, the main equations are written for this gauge as (now a dot means 



Ha 



-ul - ni - ni + + 2nin2 + 20103 + 20303 

-^[Ua(/3i,/32,/?3)], 



d 








d 






d 






1 F 

4 F20 





with Ci = 1287rGpi. 

The hamilton equation, for all Bianchi Class A cosmological models are 

III 
112 

ns 
o^ 



k = ^[-01 + 02 + 03], 
P2 = ^[-02 + 01 + 03], 

$3 = i[-03+Oi+02], 



1. Bianchi II 



01 = 2e^'^i 

02 = 0, O2 = P2 = cte, 

03 = 0, O3 = p3 = cte, 

1 F 2 

= 4 i^"^' 

A = ^ [-O1 + P2 + P3] , 

P2 = ^ [-P2 + Hi + P3] , 

P3 = J [-P3 + Hi + P2] , 

■ 1 ^ 
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introducing ()65p into (I56p we find the differential equation for Hi as Hi = — + bHi + c where 

- I (P2 + pI + Ci) • The solut ion for Hi is 



the constants are defined as b = p2 + ps and c 



Hi = b + \/-b2 - 2cTan 
and the solutions for (3i then are 

A/3i = 



^\/-b2 -2cAt 



4 Log 



Cos QV-b2-2cAtj 



A/32 = ^psAt + ^Log 
A/33 = ^P2At + ^Log 



Cos Q\/-b2 - 2cAt^ 
Cos QV-b2-2cAtj 



(73) 

(74) 

(75) 

(76) 
(77) 



and the solution for the (t> field is similar to 



-F(</.)</.2 = <Po, 
So, the solutions in the original variables are 



^dt, 



n 

/3+ 



1 

6 
V3 

6 

1 

12 



(P2 + Ps) At + Log 



Cos (^^^-b2-2cAtj 



--psAt - Log 



Cos (^i^-b2 - 2cAt j 



(p3 - 2p2) At - 2Log 



Cos QV-b2 - 2cAt^ 



(78) 



(79) 



2. Bianchi VIh=-i 



Hi 

n2 

^3 



4g2(/3i+/32) 
4g2(/3i+/32) 



0, 

1 F 



n2 = ni + ai, 

Hs = P3 = cte, 



4 F2(/) 



- [-Hi + n2 + Ps] , 
i [-n2 + Hi + Ps] , 
I [-PS + Hi + n2] , 



2F 



(80) 
(81) 
(82) 
(83) 

(84) 
(85) 
(86) 
(87) 
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introducing (I8ip into ()56p we find the differential equation for IIi as Hi — p^Hi + A;i = where 
ki = i (p| + af — 1600 + Ci — 2aiP3) who solution become as 

ni = - [eP^^* + ki], (88) 

P3 ^ ^ 

then the solutions for /3j become 

A/3i = ^(ai+p3)At, (89) 
A/32 = ^(p3-ai)At, (90) 
A/33 = 7(ai - P3)At + -L [eP3At + ^i] , (91) 

4 2p3 

(92) 



and the solutions in the original variables are 

1 



n = [2ki + p3 (ai + p3) At + 2eP3^*] , 

12p3 

/3+ = --^ [2ki + p3 (ai - 2p3) At + 2eP^^*] . (93) 



V. QUANTUM SCHEME 



The WDW equation for these models is achived by replacing IlqM = —id. 



M in (ETl). The factor 



e may be factor ordered with ILq in many ways. Hartle and Hawking 27[ have suggested what 
might be called a semi-general factor ordering which in this case would order e"^^!!^^ as 

_e-(3-Q)na^e-Q^9n = -e'^^ + Qe-'^dn, 
6 ,s d d 6 ^2 6s,_i d 

where Q and s are any real constants that measure the ambiguity in the factor ordering in the 
variables 0, and cf). We will assume in the following this factor ordering for the Wheeler-DeWitt 
equation, which becomes 

° * - V + 7* 8? + "^sn - &) * - ci* = 0. (85) 

where □ is the three dimensional d'Lambertian in the = (r2,/3+,/3_) coordinates, with signature 
(- + +)• 
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When we introduce the Ansatz ^ = xi<P)'ipi^, /3±) in (|95p . we obtain the general set of differential 
equations (under the assumed factor ordering) for the Bianchi type IX cosmological model 



6 d^x 6s 2 



(96) 
(97) 



When we calculate the solution to equation (i97|l . we find interesting properties on this, as 

1. This equation is a master equation for the field (j) for any cosmological model, implying that 
this field (j) is an universal field as cosmic ground, having the best presence in the stiff matter 
era as an ingredient in the formation the structure galaxies and when we consider two types 
of functions, 'P{4)) = uxp"^ and F{(p) = uje^'f', we have the following exact solutions 28] 



(a) F((^) = wc/)™ 

the differential equation to solver is 



with a = The solutions depend on the value to m and s, 



(98) 



i. General solution for any m ^ —2 and s, are written in terms of ordinary and 
modify Bessel function, 



1 + 3 / 2-v/a m + 2 

Vm + 2 



with ci an integration constant, Zi, is a generic Bessel function, ly 



l+s 
m+2 



(99) 



is the 



order. When a > imply uj > Q, become the ordinary Bessel function, ( Jj^, Yy). 
If a < 0, ^ t(; < 0, ^ (4, ET,,). 
ii. m = —2 and any s, 

c\ 0^ + C2<?!)~'^ si /i > 

ci + C2Ln(^ si /i = (100) 

ci sin (uLrKp) + C2 cos {fiLncp) if fi < 



X 



l + s 
) 2 < 



where = ^ -y/Kl + s)^ — 4a| 
iii. m 



—6 and s = 1 



ci sink 



"2^ 



+ C2 COs/l 



20^ 



ci sm I ^ ] + C2 cos ( ^ 



a < w < 
a > w > 



(101) 
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(b) F{(l)) = Loe"^'^, for this case we consider the caso s = 0, 

^ + ae^^X = (102) 

i. m 7^ 

X = CZo('-^e^) (103) 
\ m J 

with C is a integration constant and Zq is the generic Bessel function to zero 
order. So, if a > then u > 0, Zq is the ordinary Bessel function (Jo,lo)- When 
a < 0, w < 0, Zo {Io,Ko). 

ii. for m = 0, 



ci sinh |a|(/'j + C2 cosh (^y |a|i?i'j ifa<0— )-a;<0 
ci sin ( \/|a|(/' ) + C2 cos ( v|a[i;^ ) ifa>0— ;'a;>0 



2. If we have the solution for the parameter s=0 for arbitrary function F(i;^), say XO: then we 
have also the solution for s=-2, as xi^ = —2) = 

To obtain the solution of the other factor of ^ we use the particular value for the constants 
Ci = /i^, and make the following Ansatz for the wave function 

i;{£^) =W{i^)e-^^^"\ (105) 

where S{i^) is known as the superpotential function, and W is the amplitude of probability to 
that employed in Bohmian formalism those found in the literature, years ago 22]. So ([96|) is 
transformed into 

dW dS r 9 1 

aw -was -2VW -vs + Q— -QW— + W {vsf-u =o, (loe) 

where □ = G^^J^, V W • V $ = G^'^f^l^, (V)^ = G'^'^ J. = f + {^ f + {^)\ 
with G^^^ = diag(— 1, 1, 1), U is the potential term of the cosmological model under consideration. 
Eq (|106p can be written as the following set of partial differential equations 

{VSf-U = 0, (107a) 

W (^S + q|^^ +2VW -V S = , (107b) 

dW 

UW + Q— = 0. (107c) 

Following reference [s^, first we shall choose to solve Eqs. (|107a|) and ()107bp . whose solutions 
at the end will have to fulfill Eq. p07cp . which play the role of a constraint equation. 
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A. Transformation of the Wheeler-DeWitt equation 



We were able to solve ()107ap . by doing the change of coordinates (f22]l and rewrite ()107ap in 
these new coordinates. With this change, the function S is obtained as follow, with the ansatz 
(UMI), 

In this section, we obtain the solutions to the equations that appear in the decomposition of the 
WDW equation, (I107ap . (1107bp and (1107cp . using the Bianchi type IX Cosmological model. So, 
the equation [V]^ = — (^)^ + (gf;]:)^ + tie written in the following way (see appendix 

section 9) 



_d_ 



d d d d d d 



. d d d 



12 



d(3i d(32 dPi a/33 9/32 9/33 



d d d d d d 



dhdp2 d(3id/33 dp2dp3 
The potencial term of the Bianchi type IX is transformed in the new variables into 

2 



(108) 



U = 12 



(109) 



Then (jl07ap for this models is rewritten in the new variables as 



^',9^1+9/32 + 9/33; 



95 95 95 95 95 95 



9/3i 9/32 9/3i9/33 d^d^z 



12 



e2^l + ^2/32 + g2/33 



4e 



2(/3i+/32) _ 4g2(/3i+/33) 



4e 



2{/32+/33) 



0. 



(110) 



Now, we can use the separation of variables method to get solutions to the last equation for the S 
function, obtaining for the Bianchi type IX model 



Six = ±(e2^^+e2/^^+e2/^^ 



(111) 



In table IV we present the corresponding superpotential function S and amplitude W for all Bianchi 
Class A models. 

With this result, and using for the solution to ()107bp in the new coordinates /3j, we have for W 
function as 



Wjx = Woe[(^+t)('^i+'^2+^3l 



(112) 



and re-introducing this result into Eq. (jl07cp we find that Q = ±6. Therefore we have two wave 
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functions 



V'ix(A) = Wjx(A)Exp 



2/32 , «2/33 



Exp 



2/32 



,2/33 



Wo, Q=-6, 
WoExp[2(/3i + /32 + /33)], Q=6 



(113) 



31[ in standard quantum cosmology in general 



similar solutions were given by Moncrief and Ryan 
relativity. In table IV we present the superpotential function S, the amplitude of probability W 
and the relations between the parameters for the corresponding Bianchi Class A models. 



Bianchi 
type 


Superpotential S 


Amplitude of probability W 


Constraint 


I 


constant 


e(§ + F + 4^)A + (3 + t-4^)/52 + (3-t)/33 


r'^ -Qr-a'^ ^ 0, 
a^ = b^ + 


II 


e2/3i 


g(a-l-f)/3i+a/32 + (a-b)/33 


144b2 - 144ab + 36 
-Q2 + 24aQ = 


VIh=-i 


2(/3i -/32)e('5i+'5=) 


ea(/3l+/32) 


= 


VIIh=o 




e(l+§)(/3i+/32+/33)+a(0i+&) 


Q2- 48a -36 = 


VIII 


^2/3i g2,32 _ g2/33 




Q = ±6 


IX 


e2,3i ^ g2/3. ^ g2/33 


Wne['^+^)(^i+'3=^+'^^)] 


Q = ±6 



Table IV. superpotential S, the amplitude of probability W and the relations between the parameters for the 

corresponding Bianchi Class A models. 



If one looks at the expressions for the functions S given in table IV, one notes that there is a 
general form to write them using the 3x3 matrix m^ that appear in the classification scheme of 
Ellis and MacCallum 



32] and Ryan and Shepley [33(, the structure constants are written in the 



form 



C]k = ejksin'^^ + ^lkaj]. 



(114) 



where aj = for the Class A models. 

If we define gi(/3i) = (e'^^ , e'^^ , e^^ ) , with /3i given in (I22p . the solution to (I107ap can be written 

as 



S(A) = ±[giMy(gj)T]. 



(115) 



20 



where M'J = m'J for the Bianchi Class A, excepting the Bianchi type VIh=-i for which we redefine 
the matrix to be consistent with (IllSp 

'^0 1 o\ 

= -P2) 10 

yo oy 

Then, for the Bianchi Class A models, the wave function ^ can be written in the general form 



^ = x(</')W(/3i)exp [±[giM'j(gj)^ 



(116) 



VI. FINAL REMARKS 



Using the analytical procedure of hamilton equation of classical mechanics, in appropriate co- 
ordinates, we found a master equation for all Bianchi Class A cosmological models, we present 
partial result in the classical regime for three models of them, but the general equation are shown 
for all them. In particular, the Bianchi type I is complete solved without using a particular gauge. 
The Bianchi type II and VIh=-i are solved introducing a particular gauge. An important results 



yields when we use the gauge N = e 



,/3i+/32+/33 



, we find that the solutions for the (j) field are inde- 



pendent of the cosmological models, and we find that the energy density associated has a scaling 
behaviors under the analysis of standard field theory to scalar fields 



23 



2j], is say, scales exactly 



as a power of the scale factor like. 



oc a 



*. More of this can be seen to references cited before. 



On the other hand, in the quantum regime, wave functions of the form ^! = We^^ are the only 
known exact solutions for the Bianchi type IX model in standard quantum cosmology. In the SB 
formalism, these solutions are modified only for the function Xi ^ = xl*/") e^^^^''^ when we 

include the particular ansatz Ci = /i^. This kind of solutions already have been found in super- 



symmetric quantum cosmology 
of the heterotic strings 



34( 1 and also for the WDW equation defined in the bosonic sector 



361 ] appears all solutions in the supersymmetric 



35l | . Recently, in the books 
scheme similar at our formalism. We have shown that they are also exact solutions to the rest 
of the Bianchi Class A models in SB quantum cosmology, under the assumed semi-general factor 
ordering ([M|) . Different procedures seem to produce this particular quantum state, where S is a 
solution to the corresponding classical Hamilton- Jacobi equation (|107ap . 
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VII. APPENDIX: ENERGY MOMENTUM TENSOR 

From Eq. ([6]) we see that the effective energy momentum tensor of the scalar field is 

this energy momentum tensor is conserved, as follows from the equation of motion for the scalar 
field 




(118) 



Now we proceed to show that the energy momentum tensor has the structure of an imperfect stiff 
fluid, 

Ta /3 = (P + P)U„U^ + Pga /3 = (2/>) [U„U;3 + ^g„ ^] (119) 

here p is the energy density, p the pressure, and Ua the velocity If we choose for the velocity the 
normalized derivative of the scalar field, assuming that it is a timelike vector, as is often the case 
in cosmology, where the scalar field is only time dependent 



= S-i/2^,„, S = -</),.</)'^ (120) 

It is evident that the energy momentum tensor of the SB theory is equivalent to a stiff fluid with 
the energy density given by 

^2 2 ■ ^ ' 

Therefore the most important contribution of the scalar field occurs during a stiff matter phase 
that is previous to the dust phase. 

VIII. APPENDIX: OPERATORS IN THE A VARIABLES 



The operators who appear in eqn (I95p are calculated in the original variables (fJ, /3_|_, /3_); how- 
ever the structure of the cosmological potential term gives us an idea to implement new variables. 
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considering the Bianchi type IX cosmological model, these one given by eqn (j22p . The main 
calculations are based in the following 



d_ 

d 

9/3+ 

d 

dpi 



d d d 

+ TTT^o + TTT^o +2 



di 



d d 

+ 



9/31 

2A 

a/33' 



92 



+ 



92 



+ 



92 



9/3i9/32 9/3i9/33 d^^dh 



9/32 
%/3 



5' o 

+ 4—77 + 2 



9/3| 



'9/3| 



92 



92 



92 



9/3i9/32 9/3i9/33 9/329/33 



9 

9A 
92 

+ 



9/32 



9 

"9^ 
9/32 " 



92 



9/3i9/32 



So, the operator (V)2, □, VS'VPF are written as 



(V)^ 



9 9 

9^^ 9^^ ' 
2 



G'^'^ = diag(-l,l,l) 




+ 



9 

9^ 



+ 



9 

9^ 



^'^ = (Q,/3+,^i) 
9 9 9 9 



9 9 9 
9^ + 9;5^ + 9;5^ 



□ 



VS- VW 



92 



, 9S 9W 



92 



9/5i2 



+ 



9/32 



+ 



9 9 
9^9^ 
92 



9/3i 9/32 
9 9 



+ 



+ 



9/32 



6 



9/3i 9/33 
92 



+ 



9/3i 9/33 
9 9 



+ 



9 9 

9;s^9;s^ 



+ 



9/32 9/33 
92 



+ 



92 



9/3i9/32 9/3i9/33 9/329/3; 



9S 9W 9S 9W 9S 9W 
+ — — + 



^9/3i9/3i 9/32 9/32 9/33 9/33 
9S 9W 9S 9W 9S 9W 

+ TTT^TTT^ + 



+ 



9S 9W 9S 9W 9S 9W 

+ TTT^TTT^ + 



9/3i 9/32 9/3i 9/33 9/32 9/33 9/32 9/3i 9/33 9/3i 9/33 9/3; 



(122) 



(123) 
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